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1. Introduction 

Throughout this work, we shall consider a real- valued Levy process ^ = (Ct, i > 
0), refering to [511201 for background. This means that the process ^ starts from 
^0 = 0, has right-continuous paths with left-limits, and if {Tt)t>o denotes the 
natural filtration generated by ^, then the increment £,t+s — is independent of 
J^t and has the same law as for every s,t >0. It is well-known (see e.g. Sato 
lOnj) that the distribution of ^ is determined by its one-dimensional marginals, 
and thus by its characteristic function which has the form 

E (exp(iA^t)) = exp(-t*(A)) , A e M , t > 0. 

The characteristic exponent \E' is given by the celebrated Levy-Khintchine 
formula 

*(A) = iaA + la^X^ + / (l - e^^" + iAxl{|,|<i})n(dx) , (1) 

where a G R, cr^ > is known as the Gaussian coefficient, and 11 is a positive 
measure on IR\{0} such that /(I A |a;p)n(da;) < oo which is called the Levy 
measure. In turn, a, and 11 are uniquely determined by 4". As an example, 
for = aBt + bt where {Bt,t > 0) is a standard Brownian motion, we have 
\1'(A) — —iXb + icT^A^, so a = — 6 and 11 = 0. In the sequel, the trivial case 
when £, = (i.e. = 0) will be implicitly excluded. 

*This is an original survey paper. 
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In this text, we will be interested in the exponential functional 




exp(-6)ds , 



t > 



(note the — sign in the exponential) and especially in its terminal value 



under conditions on ^ which ensure the finiteness of loo- Our motivation stems 
from the fact that this functional arises in a variety of settings (self-similar 
Markov processes, mathematical finance, random processes in random environ- 
ment, Brownian motion on a hyperbolic space, ...), see the forthcoming Sec- 
tions [S] and El The present paper is intended as a survey of this area; it does not 
contain new results, and only a selection of the statements are carefully proven 
(of course we shall provide precise references for the remaining ones) . 

The rest of this survey consists of five sections, which are devoted to the 
following topics: Section [3 provides necessary and sufficient conditions for the 
finiteness of lac. In Section |21 we give, under certain assumptions on the Levy 
process, explicit formulae for the entire moments of the exponential functional, 
which often determine uniquely its distribution. In the general case, this dis- 
tribution can be rather complicated; nonetheless, in Section 01 we discuss some 
instances for which it can be described explicitely, in particular when ^ is a 
multiple of the Poisson process. In some cases also, the law of loo is abso- 
lutely continuous with respect to Lebesgue measure, and its density satisfies an 
integro-differential equation. Decompositions of loo before and after a stopping 
time give rise to solutions of random afhne equations. Self-similar Markov pro- 
cesses, which form an important family of Markov processes on ]0, oo[ having 
close connections with the exponential functional, are discussed in Section 
Finally, we briefly present in Section several areas where results on the expo- 
nential functionals or self-similar Markov processes can be fruitfully applied. 

2. Finiteness of the terminal value 

As a preliminary to the study of the exponential functional, wc first investigate 
when its terminal value is finite. The following result relates the finiteness of 
loo to pathwise properties of the Levy process ^, which in turn are specified in 
terms of the one-dimensional distributions of ^ or of its characteristics. 

Theorem 1 The following assertions are equivalent: 

(i) loo < oo a.s. 

(ii) P(/oo < oo) > 0. 

(iii) ^ drifts to +oo, i.e. limt^oo = +oo a.s. 

(iv) limt_+ooi"^6 > a.s. 




(v) P(^t < 0)t-i dt <oo. 
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(vi) Either 

/ \x\Il{dx) < oo and -a+ / xll{dx) e]0, oo] , (2) 

J]-oo,-l\ J\x\>l 



or 



L 



|a;|n(c?a;) = / xll{dx) = oo and I 11 {x)d{x/ J^{x)') < oo , 

oo,-l[ "'ll,oo[ -li 



where for every a; > 

n+(x) = n(]x,c»[) , n"(x) = n(]-oo,-x[) , j+ix) = f n+(y)dy. 

Jo 

Remark. The first part of condition ^ is equivalent to requiring that the 
negative part of fi, := |fi|l{^i<o}: has a finite mean. Then — — has 
a well-defined (possibly infinite) mean 

m E(Ci) = E(Ci^) -E(Cr) = -a+ / xn{dx) , 

J\x\>l 

and the second part of (0) amounts to requiring that m > 0. 
Proof: It is plain that (iv) (i) => (ii). That (ii) (i) follows from the 0-1 
law for ^. The equivalence (iii) 4^ (v) is the criterion of Rogozin (see Theorem 
VI. 12 in whereas (iii) (iv) <^ (vi) is a version of Erickson's theorem for 
Levy processes (see [S])- Finally, it is known (see Theorem VI. 12 in [S]) that 
when (iii) fails, ^ visits the negative half-line at arbitrarily large times a.s. It is 
easily seen from an application of the strong Markov property that then 1^0 ~ oo 
a.s., which entails (ii) => (iii) and thus completes the proof. ■ 



3. Calculation of entire moments 

In this section, we are interested in the calculation of the entire moments of the 
exponential functional evaluated at some independent exponential time T, viz. 
E(/|.) for fc G Z. More precisely, positive entire moments are computed when 
the Levy process ^ is a subordinator, whereas negative entire moments are being 
dealt with when the Levy process has finite positive exponential moments. 

3. 1 . Positive moments in the case of subordinator s 

In this section, we suppose that ^ = (^t, i > 0) is a subordinator, i.e. ^ has non- 
decreasing sample paths. Equivalently, this means that the Gaussian coefficient 
cr^ is equal to and that the Levy measure does not charge ] — cxo, 0] and fulfills 
/(I A x)Ii{dx) < oo, and finally 

d := -a- / xn{dx) > 0. 

^[0,1] 
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The coefficient d above is known as the drift of the subordinator. 

In this setting, it is more convenient to use Laplace transforms rather than 
characteristic functions to specify distributions. More precisely, the law of the 
subordinator ^ is characterized via the Laplace transform 

E(exp(-g6)) - exp(-i$(<z)), t,q>0, (3) 

where $ : [0, oo[^ [0, oo[ denotes the so-called Laplace exponent (sometimes also 
referred to as the Bernstein function) of ^. In turn, $ is given by an analytic 
continuation of the characteristic exponent ^' to the upper half-plane, namely 

^{q) = l'(iq), q>0. 

The Levy-Khintchine formula now reads 

^{q)=dq+[ (l-e-9^)n(dx), q>0. 
J]0,oo[ 

We refer to [7j for background on subordinators. 

The first part of the following theorem has been established by Carmona et 
al. 1201 J the second is an immediate extension of 15 . 

Theorem 2 Suppose that ^ is a subordinator with Laplace exponent $. Let T 
denote a random time which has an exponential distribution with parameter 
9 (for q = 0, we have T = oo) and is independent of ^. 

(i) The positive entire moments of the exponential functional evaluated at time 
T are given in terms of $ by the identity 

As a consequence, E(exp(a/T)) < oo whenever a < ^(oo), and the distribution 
of It is determined by 

(ii) For each q > 0, there exists a unique probability measure pq on [0, oo[ which 
is determined by its entire moments 

x'^Pgidx) ^ iq + $(1)) ■ • ■ (g + $(fc)) fork = 1,2,... 

[0,oo[ 

In particular, if Rq is a random variable with law pq that is independent of ^, 
then we have the identity in distribution 

IrRq = e , (5) 
where e denotes a standard exponential variable. 

Of course, one can rephrase formula ^ in terms of Laplace transform for the 
moments of the exponential functional evaluated at a time i, specifically 

/ £(1^)0-9*^*=- TTT—, ^TTTT, for fc = 1, 2, . . . 

7o ^ ' 5(g + $(l))---((7 + $(fc)) 
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which provides a formal characterization of the distribution of It. However in 
practice, an exphcit inversion of the Laplace transform does not seem to be an 
easy matter. 

Proof: (i) We first establish the formula Q in the case q — Q, that is for 
T ~ oo a.s. Set 

/oo 
exp{-£_s)ds 

for every t > 0. On the one hand, for every positive real number r > 0, we have 
the identity 

Jo-JI = r I eM-is)Jr'ds. (6) 
Jo 

On the other hand, we may express Jg in the form Js = exp(— ^s)/^, where 

/"OC 

I'^ = / exp{-Odu and C - 6+« - 6 • (7) 

From the independence and stationarity of the increments of the Levy process, 
we see that I'^ has the same law as loo and is independent of ^g. Plugging this 
in © and taking expectations, we get using Q that 

E(/^)(l-cxpH$(r))) ^ rj\M~s^r))Eir^')ds 



Finally 



and since E(7^) = 1, we get the formula of the statement for q = by iteration, 
taking r = fc G N. 

Now we treat the general case when the parameter q of the exponential time 
is finite. Pick a > arbitrarily large, and let N — {Nt,t > 0) be an independent 
Poisson process with parameter q. So if T denotes the first jump time of the 
Poisson process N, then T has an exponential distribution with parameter q and 
is independent of ^. Write ^1°''' = + aNt, so is again a subordinator with 
Laplace exponent ^^"-^{p) = + ^(l - e~'^P), p >0. Denote also by /'"^ the 
exponential functional of Plainly, J^'' decreases to It as a increases to oo, 
and since lima^oo <l'^°-'(p) = ^(p) + q for every p > 0, we deduce by dominated 
convergence that 

E(/^) = lim £((4"^)'=^) 



= lim 



T 

kl 



a^oo $(1) (1) • • • (fc) 

fc! 

(g + $(!))• ••(g + $(fc)) 
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The rest of the proof is immediate. 

(ii) We refer to for the proof of the existence of p; see also Berg and 
Duran ^ for an analytical proof. The second part of the statement now follows 
from (i), since k\ is the fc-th moment of the standard exponential law and the 
latter is determined by its entire moments. ■ 

In Section 14.11 we shall present a couple of examples in which the moment 
problem Q can be solved explicitly, yielding expressions for the law of the 
terminal value of exponential functional /. Theorem |2Iii) also invites to look 
for factorizations of an exponential variable in two independent factors, one of 
which may be obtained (in distribution) as the exponential functional loo- We 
refer to |E| for several classical such factorizations which fit in this setting. 

We conclude this section by mentioning that moments calculations as in The- 
orem|2Ii) can be extended to the situation when ^ is an increasing additive pro- 
cess, i.e. its increments are independent but not necessarily homogeneous. This 
situation arises naturally in particular in the representation of self-decomposable 
laws on M_|_ (see Sato (55] and Jeanblanc, Pitman and Yor |41j)- It also occurs in 
in non-parametric Bayesian statistics, where the mean of a random distribution 
chosen from a neutral-to-the-right prior can be represented as the exponential 
functional of an increasing additive process. See |28| and the references therein 
for much more on this topic. 

3.2. Negative moments 

Here, we shall suppose that the Levy process ^ possesses exponential moments 
of all positive orders, i.e. E(exp(g^t)) < oo for every t,q > 0. Recall (e.g. from 
Theorem 25.3 in Sato (6|)|) that this is equivalent to assuming that the Levy 
measure fulfills 

/ e'^^n(da;) < oo for every q > 0; 

J[l,oo[ 

this condition holds for instance when the jumps of ^ are bounded from above 
by some fixed number, and in particular when ^ has only negative jumps. In 
this situation, the characteristic exponent 5* has an analytic extension to the 
half-plane with negative imaginary part, and one has 

IE(e««') = exp(tV(g)) < oo, t,q>0, (8) 

where 

^{q) := -*(-ig), q>0. 
In this direction, note that the assumption that ^ drifts to +oo, is equivalent to 

m = E(ei) - ^-'(0+) e]0,oo[, (9) 

see the remark after Theorem ^ The following result appears in |15|. 
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Theorem 3 Assume that ^ and ^ hold. Then, for every integer k > 1 and 
for every < t < oo, E(/j^'^) < oo; furthermore we have 



, _ V.(i)...V.(fc-i) 

^^""^""^ (fc-1)! 



with the convention that the right-hand side equals m for k ^ I. 

If moreover ^ has no positive jumps, then l/Ioo admits some exponential 
moments. Hence, the distribution of loo is determined by its negative entire 
moments. 

Proof: To start with, let us first check the finiteness of the negative moments 
of It for t > 0, and a fortiori for loc- In this direction, we observe from JH)) that 
for every a, g > 0, 

P(a > a) < exp(^(<z) - qa) . 
Next, we have for all sufficiently large a that 

P( inf 6 < -a) < 1/2, 

0<t<l 

and it follows from an application of the strong Markov property at the first 
passage time of ^ above the level 2a that 

P( sup ft > 2a) < 2P(fi > a) < 2exp{^jj{q) - qa) . 

0<t<l 

This entails the bound 

Wi < < 1P( sup ft > 21ogl/£) < 2exp(V'((?))£« 

0<t<l 

for small enough £ > 0, and thus ensures that E(/j~'^) < oo for every fc > 0. 

Now the proof is similar to that of Theorem I^Ji) . More precisely, recall the 
notation 



/oc 



for every t > and observe that for every real number r > 0, there is the 
identity 

Jr - Jo"' - r f cxp(-f,) J7('-+i)ds . (10) 
As in (0, we write for every fixed s 

Js = exp{-£,s)I'oo , 

where has the same law as /qo and is independent of f^. Plugging this in 
l(TU|l and taking expectations, we get from ||5J that 

E (J--^) (exp(i^(r)) - 1) = r T exp(si/;(r))E(/-(-+i)) , 
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and finally 

(11) 

Because we know that the negative moments of lac are finite, we may let r — > 0+ 
and apply the theorem of dominated convergence. We get 

IE(l//oc) ^ lim = ^'{0+) = m. 

r^0+ r 

The formula of Theorem|21now follows by induction. 

Finally, let us check that the distribution of 1 /loo is determined by its (pos- 
itive) entire moments when ^ has no positive jumps. In that case, the Levy 
measure IT does not charge [0, oo[, and the Levy-Khintchine formula ^ readily 
yields the bound ipil) = 0{q^) as g ^ cxd. As a consequence, there is some finite 
constant c > such that 

£(/-'=) = m ^^^^^'^'^l^'^^ < c^kl for every fc = 1, . . . . 
This entails that 



E 



exp 



1 



< oo , 



hence the distribution of \/Ioo is characterized by its entire moments. 



We make the conjecture that the positive entire moments of l/Ioo do not 
determine its distribution when the Levy process ^ has positive jumps (it has 
been shown in |T3 that this is indeed the case - discussed in the next subsection 
- when ^ is a Poisson process, and more precisely, that the negative entire 
moments of the exponential functional of a Poisson process are closely related 
to those of the log-normal distribution). We refer to |fi2l lf)4j and the references 
therein for more about indeterminacy for the moment problem. 

4. On the distribution of the exponential functional 
4.1. Some important special cases 

In this section, we briefiy present some explicit results on the distribution of the 
exponential functional in three important special cases; further examples can 
be found e.g. in 20, 21 . 

• Standard Poisson process: Here, we consider a standard Poisson process 
{Nt,t > 0) and define for q g]0, 1[ its exponential functional by 



Jo 
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In other words, — loo for ^ — — {log q)N. Note that we may also express 
l!^^ in the form 

oo 

= E 9"^" > (12) 

where = Tn+i — T„, n — 0, 1, . . . denote the waiting times between the 
successive jump times T„ = mf{t : Nt = n} of {Nt,t > 0). In other words, 
(e„,n G N) is a sequence of i.i.d. exponential variables with parameter 1. 

In order to specify the distribution of the exponential functional, it is quite 
convenient to use the so-called g-calculus (see, e.g. 03, •.•) which is asso- 
ciated with the basic hypergeometric series of Euler, Gauss, ... More precisely, 
let us introduce some standard notation in this setting: 

71 — 1 OG 

(a;g)« = ]^(1 - ag^) , (a; g)oo = ^(1 - ag^) , 
j=o J=0 



and 



We may now state the main formulas for the distribution of the exponential 
functional of the Poisson process, which are excerpts from ^U] (see also |25|1: 
The Laplace transform of 7^'' is given by 



E (exp(A4^))) = j—^ (A < 1) , (13) 
its Mellin transform by 



and its density, which we denote as {i^'^^{x), x > 0), by 

z(')(x) = E (-^Z*^") / , , , • (15) 

^ {q;(i)oo[q;q)n 

More precisely, follows immediately from H12|) . whereas (|14|l can be de- 
rived from the factorization (|SJ of the exponential law. Finally, \ib\ is obtained 
by inversion of the Mellin transform; see 10 for details. 

• Brownian motion Vifith drift: Take = 2{Bt + bt), where is a standard 
Brownian motion and & > 0. Then m = 26, ipil) — 2q(g -I- b), and we get for 
every integer fc > 1 

mm_l^ = 2^-6(^+1). ..(6 + fc-l) 

gfc r(fc + b) 

m ■ 
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The right-hand side can be identified as the k-th moment of 27b, where 7b is 
a gamma variable with index b. Hence Theorem |31 enables us to recover the 
identity in distribution 

/ exp {-2(B, + bs)} ds = (16) 

Jo ^Ib 

which has been established by Dufresne I^HI (see also Proposition 3 in PoUak- 
Siegmund 1^7 , Example 3.3 on page 309 in Urbanik [HS], Yor ...). A further 
discussion of (|16|1 in relation with DNA is made in 0^ and |48j . 

• Exponential Levy measure: Suppose ^ has bounded variation with drift 
coefficient +1 and Levy measure 

T\{dx) = {a + b- l)5e^^dx , x<0, 

with 0<a<l<a + b. Then, ^(g) = q{q + 1 - a) / (b + q), m = (1 - a)/b, and 
we get for every integer A: > 1 

(l-a)---(fc-a) 
(fc-1)! 6---(fc + 6-l) 

r(A: + 1 - a)r(6) 

r(i -a)r(fc + 6) ■ 

On the right-hand side, we recognize the fc-th moment of a beta variable with 
parameter {1 — a,a + b — 1). We thus get the identity in law 

c 1 

T — 



Pl-a,a+b-l 

which has been discovered by Gjessing and Paulsen [5^ . 
4-2. An integro- differential equation 

In this section, we present an equation satisfied by the density of the law of the 
exponential functional. We shall assume here that the Levy process is given in 
the form 

= ct + aBt + Tt+ - TfT , 

where (_Bj, i > 0) is a standard Brownian motion, c a real number, two subor- 
dinators with no drift, such that B, t+ and are independent. Using the same 
notation as in Thcorcm^vi), we shall further suppose that e^II~(x)da: < 00 
and that 

/•oo /"OO /"OC 

/ Il'^{x)dx < 00 , c+ n+(a;)dx- / Il^{x)dx < 00, 
Jo Jo Jo 

(which is equivalent to assuming that ^1 is integrable and has a positive mean, 
cf. the remark after Theorem^. Carmona et al. 21 have proved the following 
result. 
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Theorem 4 Let the above assumptions prevail. Then the law of the terminal 
value loo of the exponential functional is absolutely continuous with respect to 
Lebesgue measure and has a C°° density k which solves the following eguation: 



a ^ 2 



^^^ix^kix))+[y-~cjx + ljkix) 

IV+ {\og{u/x))k{u)du- I IT {\og{x/u))k{u)du. 
Jo 

We refer to Carmona et al. |23 for a proof and for some applications of this 
equation. 



4.3. Random affine equation 

The strong Markov property of Levy processes has an interesting consequence 
for the exponential functional. Indeed, consider some finite stopping time T in 
the natural filtration {!Ft)t>o of ^, so that the process := £_T+t ~ Ct, i > 0, is 
independent of J-t and has the same law as ^. Now writing 

/•oo poo 

/ e^p{-Cs)ds = exp(-CT) / ejip{-Qdr + / exp(-^s)ds 
Jo Jo Jo 

yields a so-called random afHne equation 

loo = AIoo+B (17) 

where, on the right-hand side, the variables loo and {A, B) are independent. 
Specifically, the law of the latter pair is given by : 



iA,B) ^ l^exp(-eT),^ exp(-6 



)ds 



The random affine equation has been investigated in depth by Kesten j see 
also Vervaat [HS! and Goldie 

An interesting application of the results of Kesten and Goldie in the special 
case when we take T = 1 is the following estimation of the tail distribution of 
the exponential functional under a condition of Cramer's type. 

Corollary 5 Suppose there exists 9 > such that E(exp(— 0^f)) 1 for some 
(and then all) t > 0. Assume further that E(^t exp(— < 00 and that ^ is 
not arithmetic (i.e. there is no real number r > such that P{^t G ''^) = 1 for 
all t >0). Then there exists some constant c > .such that 

P(/oo >t) Ct-" , t ^ 00 . 

This interesting observation has been made by Mejane [SSI and Rivero |nHl- It 
constrasts with the situation when ^ is a subordinator (then, obviously, the con- 
dition of Corollary never holds), as we know that then loo possesses some finite 
exponential moments. Finally, we refer to the recent work |52| for extensions of 
Corollary |5l to situations where Cramer's condition is not fulfilled. 
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5. Self-similar Markov processes 

Motivated by limit theorems for Markov processes, Lamperti 07j considered 
families of probability measures {Px,x > 0) on Skorohod's space D of cadlag 
paths Lu : M+ M+, under which the coordinate process X.{lu) = w(-) becomes 
Markovian and fulfills the scaling property. That is for some real number a, 
called the index of self-similarity, there is the identity 

the law of {kX/.-cj., t > 0) under P^; is T^kx (18) 

where fc > and a; > are arbitrary. These laws are called self-similar (or 
semi-stable) Markov distributions. 

Lamperti has shown that self-similar Markov processes can be connected to 
Levy processes by a bijective correspondence that we shall present in the next 
section. Lamperti's correspondence also points at the key role of the exponential 
functional in this framework, and we shall see in particular in Section 15.31 that 
the law of /qo is fundamental to understand the asymptotic behavior of self- 
similar Markov processes. 



5.1. Lamperti's transformation 

Recall that ^ = {£,t,t > 0) denotes a real- valued Levy process which starts 
from 0, and that the notation P, E is used for probabilities and expectations 
related to ^. In the sequel, it will be convenient to use the notation 

Jo 

for every 6 G R (in other words, 7*^^-' is the exponential functional associated to 
the Levy process b^). 

For each a G R, we define a time-change r^"-* = (T("\t > 0) as the inverse 
of the exponential functional /^~"\ i.e. 

T^"^ := inf |u > : ^ exp(a^,)ds > ^| . 

More precisely, the time-change is always finite when the Levy process oscillates, 
or when ^ drifts to -|-oo and a > 0, or when drifts to ~oo and a < 0. In these 
cases, there is the identity 

I^"' = / exp(aC.)ds = t. (19) 
Jo 

When ^ drifts to -l-oo and a < 0, or when ^ drifts to — oo and a > 0, (|19|l holds 
if and only if t < /^"\ and r" = oo otherwise. 
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Next, for an arbitrary a; > 0, denote by P^; the distribution on Skorokhod's 
space D of cadlag paths of the process 

Xt := a;exp | , t>0. 

It is this transformation from Levy processes which we call Lamperti's transfor- 
mation. In the case when the Levy process drifts to +00 and a < 0, we agree 
that the quantity above equals 00 for all t > x'^lio Similarly, we agree that 
Xt ~ Q when the Levy process drifts to —00, a < and t > In other 

words, and 00 are the two absorbing boundaries; observe that these boundaries 
can only be reached continuously. 

Then, writing [Qt)^^^ for the natural filtration of the coordinate process X 
on D, we have that for every s,t > and measurable function / : R+ 

{f{Xt+s) I Gt) - EyifiXs)) for y^Xt, 

and the scaling property (|18|l holds. Conversely, any self-similar Markov process 
which either never reaches the boundary states and 00, or reaches them con- 
tinuously (in other words, there is no killing inside ]0, oo[) can be constructed 
by Lamperti's transformation from a Levy process. 

We further point that for every (3^0, the process X raised to the power f3, 
X^, is self-similar with index a/P; more precisely it can be constructed from 
Lamperti's transformation from the Levy process (3£,. This enables to focus on 
the case when the index of self-similarity is 1 . 

That exponential functionals play an important role in the study of self- 
similar Markov processes is plain in cases when absorption can occur, since the 
distribution of the absorption time under P^; is given by the identity 

C = 2:"4-") . 

We shall also see in the next sections that, at least in some special cases, there are 
explicit formulas for the moments of Xt which bear similarities with those of the 
exponential functional, and that the terminal value of the exponential functional 
/^"^ is crucial to determine the asymptotic behavior of the self-similar Markov 
process. 

5.2. Moments at a fixed time 

Throughout this section, we shall assume that a, the index of self-similarity of 
X, is positive. 

First, we consider the case when the Levy process ^ is a subordinator; recall 
that $ denotes its Laplace exponent, cf. formula Q. One gets 

00 / _p,\ji n — 1 

(Xr) = J2 n + fc") ' (20) 

n=0 ' k=0 
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where in the right-hand side, we agree that the product equals 1 for n — 0. See 
Equations (2) and (3) in |13| where this formula is proven for x = a — 1 (the 
general case follows easily) . It may be interesting to check that (|20l) agrees with 
Kolmogorov's backward equation. Specifically, if we write G for the infinitesimal 
generator of X, and fp-.x^x^ for the power function with exponent p £ 
M, then standard transformations for Markov processes (time substitution and 
image by a bijection) entail that 

It is then immediate to see that if we write ft (x) for the quantity in the right- 
hand side of jSnj, then dtft{x) = Gft{x) as required. 

We also point out that in the special case x = a = p ~ 1, formula 
shows that the function t — > Ei (l/Xt) coincides with the Laplace transform 
of the probability measure po which is related to the exponential functional in 
Theorem |2i;ii). 

Now we turn our attention the case when the Levy process has finite exponen- 
tial moments of any positive order, and which does not drift to — oo. Specifically, 
we assume that (jHJ holds and that m = E(^i) > 0. Then, following one 
can establish using Kolmogorov's backward equation that the entire moments 
of X" are given for fc = 1, ... by 

(Xr^) = x"^ + ^ ^(afc)---VXa(fc-^+l)) ^,(,,,)^, ^21) 
e=i 

Furthermore, if ^ has no positive jumps, then (|21|l determines the semigroup of 
the self-similar Markov process. 

5.3. Asymptotic behavior for positive indices of self- similarity 

Throughout this section, we shall suppose that the index of self-similarity is 
positive : a > 0, and that the Levy process ^ drifts to -l-oo. In this situation, 
recall that the self-similar Markov process tends to oo without reaching oo, and 
it is therefore interesting to obtain sharper information about the asymptotic 
behavior oi Xt as t — > oo. The following result was proved in (see also 
Bertoin and Caballero in the special case when the Levy process ^ is a 
subordinator). 

Theorem 6 Let X denote a self-similar Markov process with index a > as- 
sociated by Lamperti's transformation to the Levy process ^. Suppose that ^ is 
not arithmetic (i.e. does not live on a discrete subgroup rZ for some r > 0), 
and has a positive finite mean, i.e. G L^{¥) and lE(Ci) to > 0. Then for 
every x > 0, under P^, t~^^°'Xt converges in distribution as t ^ oo to some 
probability measure fia on R+ specified by 
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where f : R+ denotes a generic measurable function. 

Example. Suppose a — 1 and = 2{Bt + bt) for some constant 6 > 0, i.e. ^ is 
twice a standard Brownian motion with drift b. Then it is easy to check, e.g. by 
stochastic calculus, that X solves the stochastic differential equation 

dXt = 2^/x'td[3t + 2(6 + l)dt 

where (/3t , i > 0) is a standard Brownian motion. Thus X is the square of a Bessel 
process with dimension 2(6 + 1). On the other hand, we know from Dufresne's 

identity Hl()|) that /<^'' = loo = l/27b where jb denotes a gamma variable with 
parameter 6. Hence Theorem yields that the limiting distribution yUi is given 

by 

/ = 6-iE(7f,/(27f,)) = E(/(27,+i)), 

i.e. t-'^Xf converges in law as t — > oo towards twice a gamma variable with 
parameter 6 + 1 . 

Theorem El is the key to understand the entrance distribution of the self- 
similar Markov process from 0+ (i.e. the limit of the laws Px as x ^ 0+); see 
[Tl) for details and also Caballero and Chaumont for a recent development. 



6. Some applications 

Besides its role in the study of self-similar Markov processes, exponential func- 
tionals arise in a number of domains in probability theory and related topics, 
among which (the following list is by no mean exhaustive): 



6.1. Brownian motion on hyperbolic spaces 

The random variable 



exp(B('^))dW-i^) , (22) 

where Ws'^^ = Ws — fJ^s and si"'* = Bs — vs are two independent Brownian mo- 
tions with respective drifts — /i G M and —v < 0, appear both in risk theory (see 
e.g. Paulsen [21], Example 3.1) as well as in connection with invariant diffusions 
on the hyperbolic half-plane EI (see Bougerol ^Hj: Ikeda and Matsumoto j4()| . 
Baldi et al. ...). Hence, it is no surprise that the distribution of (|22|l has 
been much studied; it has a density given by 

exp(— 2/i arctanx) 

f{x) = 



(l + a;2)'^+i/2 



which belongs to the type IV family of Pearson distributions. 
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Concerning diffusions on H, consider the infinitesimal generator 

L(-) = ^A-,.|--(.-l/2),| 

{{x,y) denote the Cartesian coordinates on H). The differential operator L'^t^''^^ 
is invariant under the orientation-preserving isometrics of H that fix the point at 
infinity, that is under the real affine transformations z az + b with a > and 
6 e K. The diffusion process associated to L^^-") corresponds to the stochastic 
differential equation 

dXt = YtdWt - fiYtdt 
dYt^YtdBt-{v-l/2)Ytdt, ^^•'> 

where, as before, W and B are independent one-dimensional Brownian motions. 
The solution of (|23|l with starting point iy = (0, y) is : 

{ X, = jy/Jexp(Bi'^))dW^i^) 
\ Yt=yeMB^t'^)- 

In close connection with this diffusion, one may derive a way to obtain the 
density f{x) exhibited above. This may follow from the following three simple 
observations: 

(a) For a fixed real number x, consider the two processes 

Xl^--"^ = exp{B[-^)x + f exp(B(^))dW-i^) , 

JQ 

and ^ 

X'l'''"^ = eMB["^) + exp(-B('^^)dW^i^) 

By time-reversal, these two processes have the same one-dimensional distribu- 
tions, although they do not have the same law. 

(b) The process {x[^'^\ i > 0) is a diffusion with generator 



d 

dx^ v/- ■ \ ' ' ' dx 



^=^^33 -(M+(^-l/2)a;)- 



(c) The distribution at time t of this diffusion process converges to the in- 
variant distribution with density f{x). Indeed, from (a), x[^''^^ has the same 
law as x['^''^\ and the former converges in law as i — > cxo to, say, X^''^\ The 
limit distribution is invariant, that is annihilated by A/*, the adjoint of Af , and 
hence it coincides with /. 



J. Bertoin and M. Yor / Exponential functionals 



207 



6.2. Asymptotic studies of Brownian diffusions in random 
environments 

Exponential functionals occur very naturally in the study of some models of 
Random Walks in Random Environment, which we now describe informally. 
Associated with a stochastic process V = {V{x),x G M) such that V{0) = 0, a 
diffusion Xy — > 0) in the random potential V is, loosely speaking, a 

solution to the stochastic differential equation 

dXv{t) = d(3t - ^V'{Xv{t))dt , Xv{0) ^ , 

where {Pt,t > 0) is a standard Brownian motion independent of V. More rig- 
orously, the process Xy should be considered as a diffusion whose conditional 
generator, given V, is: 

iexp(m);i(.-"'"^). 

It is now clear, from Feller's construction of such diffusions, that the potential 
V does not need to be differentiable. 

Brox ^] and Kawazu and Tanaka |42[ 021 (see also e.g. Hu et al. 39 and 
Comtet et al. '23') have studied this random diffusion in the case when y is a 
Brownian motion with negative drift : V{x) = Bx — kx, for fc > 0. Thus, it is no 
surprise that the knowledge about the exponential functionals /q dye'xjp{By — 
ky), cc G K, plays an essential role in this domain; see Shi '61' for a nice picture 
of this area. 

6. 3. Mathematical finance 

In mathematical finance, the processes of reference were chosen originally to be 
the exponential of Brownian motions with drift (i.e. the famous Black & Scholes 
framework, which has been widely adopted between 1973 and 1990, say), and 
later generalized to become the exponential of Levy processes (Ct , * > 0) . In 
particular, the computation of the price of Asian options (see, e.g., papers in 
1201) is equivalent to the knowledge of the law of /p exp(— ^5)^5, for fixed t\ 
an important reduction of the problem consists in replacing t by an exponential 
time 0, which is independent of ^ (in other words, one kills the Levy process with 
some constant rate). Furthermore, some identities in law due to Dufresne |27| . 
for exponential functionals of Brownian motions with opposite drifts, involving 
both the case t < 00 and t — 00, led to some interesting extensions of Pitman's 
(2A/ — X) theorem for exponential hmctionals (see, e.g. |49l 1501 

6.4. Self-similar fragmentations 

Fragmentations are meant to describe the evolution of an object with unit mass 
which falls down randomly into pieces as time passes. They arise in a variety 
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of models; to give just a few examples, let us simply mention the studies of 
stellar fragments and meteoroids in astrophysics, fractures and earthquakes in 
geophysics, breaking of crystals in crystallography, degradation of large polymer 
chains in chemistry, fission of atoms in nuclear physics, fragmentation of a hard 
drive in computer science, crushing in the mining industry, ... 

The state of the system at some given time consists in the sequence of the 
masses of the fragments. Suppose that its evolution is Markovian and obeys the 
following rule. There is a parameter a G M, called the index of self-similarity, 
such that given that the system at time t > consists in the ranked sequence of 
masses mi > m2 > • . . > 0, the system at time t + r is obtained by dislocating 
every mass rm independently of the other fragments to obtain a family of sub- 
masses, say {mij,j e N), where the sequence of the ratios {mij/mi,j G N) 
has the same distribution as the sequence resulting from a single unit mass 
fragmented up to time mfr. Such a random process will be referred to as a 
self-similar fragmentation with index a. 

For the sake of simplicity, we shall focus here on the following elementary 
model. A fragment with mass m > remains unchanged during an exponential 
time with parameter m". Then it splits and gives rise to a family of smaller 
masses, msi,ms2, . . . where, (si, S2, . . .) is a random sequence of nonnegative 
integers with sum si + S2 + . . . = 1, which has a fixed distribution i^. One calls 
v the dislocation law of the fragmentation process. 

In this setting, self-similar Markov processes arise as follows. One taggs a 
point at random in the object according to the mass measure and independently 
of the fragmentation process, and one considers at each time t the mass Xt of 
the fragment containing the tagged point. Then {l/xt,t > 0) is a self-similar 
Markov process with index a, and the Levy measure 11 of the Levy process ^ 
which is associated to 1/x by Lamperti's transformation (cf. Section l5!T|l . can 
be expressed explicitly in terms of the dislocation law v as 



Statistical properties of the tagged fragment can then be derived e.g. from re- 
sults presented in Section ??, and then be shifted to the fragmentation process. 
In particular, this approach yields important limit theorems for the empirical 
distribution of the fragments of self-similar fragmentations with a positive index 
of self-similarity; see 1121 IT71 On the other hand, self-similar fragmenta- 
tions with a negative index of self-similarity are dissipative, in the sense that 
the total mass of the fragments decreases as time passes. In this direction, the 
terminal value of the exponential functional gives the absorption time of the 
self-similar tagged process, i.e. the time at which the mass of the tagged frag- 
ment vanishes. Distributional properties of the exponential functional provide a 
crucial tool for the analysis of the loss of mass; see Haas |SZ1 EHI ■ 

In a related area, we mention that the terminal value of the exponential 
functional of subordinators also arises in the study of regenerative compositions, 
see the recent work of Gnedin, Pitman and Yor [331 1^ ■ 




X €]0, Qo[ . 



J. Bertoin and M. Yor / Exponential functionals 



209 



References 

[1] P. Baldi, E.G. Tarabusi, A. Figa-Talamanca, and M. Yor (2001). Non- 
symmetric hitting distributions on the hyperbohc half-plane and subor- 
dinated perpetuities. Rev. Mat. Iberoamericana 17, 587-605. MRIOOOSQU 
[2] Ch. Berg (2004). On powers of Stieltjes moment sequence I. J. Theoret. 

Probab. (to appear). 
[3] Ch. Berg (2004). On powers of Stieltjes moment sequence II. Preprint. 
[4] Ch. Berg and A. J. Duran (2004). A transformation from Hausdorff to 

Stieltjes moment sequence. Ark. Math. 42, 239-257. 
[5] J. Bertoin (1996). Levy Processes. Cambridge University Press. 
[6] J. Bertoin (1997). Regularity of the half-line for Levy processes, Bull. Sc. 

Math. 121, 345-354. MR1465812 
[7] J. Bertoin (1999). Subordinators: Examples and Applications. Ecole d'ete 
de Probabilites de St-Flour XXVII, Lect. Notes in Maths 1717, Springer, 
pp. 1-91. MR1746300 
[8] J. Bertoin (2002). Self-similar fragmentations. Ann. Inst. Henri Poincare 

Ann. Inst. Henri Poincare 38, 319-340. 
[9] J. Bertoin (2003). The asymptotic behavior of fragmentation processes. /. 
Europ. Math. Soc. 5, 305-416. MR2017852 
[10] J. Bertoin, Ph. Biane, and M. Yor (2004). Poissonian exponential function- 
als, 5-series, g-integrals, and the moment problem for log-normal distribu- 
tions. Proceedings Stochastic Analysis, Ascona, Progress in Probability vol. 
58, Birkhauser, pp. 45-56. 
[11] J. Bertoin and M.-E. Caballero (2002). Entrance from 0+ for increasing 

semi-stable Markov processes. Bernoulli 8, 195-205. 
[12] J. Bertoin and A. Gnedin (2004). Asymptotic laws for nonconserva- 
tive self-similar fragmentations. Elect. J. Probab. 9, 575-593. Available at 
http : //www. math. Washington . edu/ ejpecp/viewarticle .php?id=1463 
&layout=abstract. MR2080610 
[13] J. Bertoin and M. Yor (2001). On subordinators, self-similar 
Markov processes, and some factorizations of the exponen- 
tial law. Elect. Commun. Probab. 6, 95-106. Available at 
http : //www.math. Washington. edu/ ejpecp/ECP/viewarticle .php?id= 
1622&layout=abstract. MR1871698 
[14] J. Bertoin and M. Yor (2002). Entrance laws of self-similar Markov pro- 
cesses and exponential functionals of Levy processes. Potential Analysis 17, 
389-400. MR1918243 
[15] J. Bertoin and M. Yor (2002). On the entire moments of self-similar Markov 
processes and exponential functionals of certain Levy processes. Ann. Fac. 
Sci. Toulouse Serie 6, vol. XI, 33-45. 
[16] Ph. Bougerol (1983). Exemples de theoremes locaux sur les groupes 

resolubles. Ann. Inst. H. Poincare 19, 369-391. 
[17] M. D. Brennan and R. Durrett (1987). Splitting intervals. II. Limit laws 
for lengths. Probab. Theory Related Fields 75, 109-127. MR879556, 



J. Bertoin and M. Yor / Exponential functionals 



210 



[18] T. Brox (1986). A one-dimensional diffusion process in a random medium. 
Ann. Probab. 14, 1206-1218. 

[19] M.-E. Caballero and L. Chaumont (2004). Overshoots of Levy processes 
and weak convergence of positive self-similar Markov processes. Preprint. 

[20] P. Carmona, F. Petit and M. Yor (1994). Sur les fonctionnelles exponen- 
tielles de certains processus de Levy. Stochastics and Stochastics Reports 
47, 71-101. English translation in [70) . 

[21] P. Carmona, F. Petit and M. Yor (1997). On the distribution and asymp- 
totic results for exponential functionals of Levy processes. In: M. Yor (edi- 
tor) Exponential functionals and principal values related to Brownian mo- 
tion, pp. 73-121. Biblioteca de la Revista Matematica Iberoamericana. 

[22] P. Carmona, F. Petit and M. Yor (2001). Exponential functionals of Levy 
processes. O. Barndorff- Nielsen, T. Mikosch and S. Resnick (editors). Levy 
processes: theory and applications, pp. 41-55, Birkhauser. MR1833689 

[23] A. Comtet, C. Monthus, and M. Yor (1998). Exponential functionals of 
Brownian motion and disordered systems. J. Appl. Probab. 35, 255-271. 

[24] R. Cowan and S. N. Chiu (1994). A stochastic model of fragment formation 
when DNA replicates. J. Appl. Probab. 31, 301-308 MR1274788 

[25] V. Dumas, F. Guillemin and Ph. Robert (2002). A Markovian analysis 
of additive-increase, multiplicative-decrease (AIMD) algorithms. Adv. in 
Appl. Probab. 34, 85-111. 

[26] D. Dufresne (1990). The distribution of a perpetuity, with applications to 
risk theory and pension funding. Scand. Actuar. J. 1-2, 39-79. 

[27] D. Dufresne (2001). An affine property of the reciprocal Asian option pro- 
cess. Osaka J. Math. 38, 379-381. 

[28] I. Epifani, A. Lijoi and I. Prnster (2003). Exponential functionals and means 
of neutral-to-the-right priors Biometrika 90, 791-808. 

[29] W. E. Feller (1971). An introduction to probability theory and its applica- 
tions, 2nd edn, vol. 2. Wiley, New York. 

[30] A. F. Filippov (1961). On the distribution of the sizes of particles which 
undergo sphtting. Th. Probab. Appl. 6, 275-293. 

[31] G. Gasper and M. Rahman (1990). Basic Hypergeometric Series. Cam- 
bridge University Press, Cambridge. 

[32] H. K. Gjessing and J. Paulsen (1997). Present value distributions with 
application to ruin theory and stochastic equations. Stochastic Process. 
Appl. 71, 123-144. 

[33] A. Gnedin and J. Pitman (2005). Regenerative composition structures. 
Ann. Probab. 33, 445-479. 

[34] A. Gnedin, J. Pitman and M. Yor (2004). Asymptotic laws for com- 
positions derived from transformed subordinators. Preprint available via 
http : //front . math . ucdavis . edu/math . PR/0403438. 

[35] C. M. Goldie (1991). Implicit renewal theory and tails of solutions of ran- 
dom equations. Ann. Appl. Probab. 1, 126-166. 

[36] F. Guillemin, Ph. Robert, and B. Zwart (2004). AIMD algorithms and 
exponential functionals. Ann. Appl. Probab. 14, 90-117. 

[37] B. Haas (2003). Loss of mass in deterministic and random fragmentations. 



J. Bertoin and M. Yor / Exponential functionals 



211 



Stochastic Process. Appl. 106, 245-277. 

[38] B. Haas (2004). Regularity of formation of dust in self-similar fragmenta- 
tions. Ann. Inst. Henri Poincare 40, 411-438. 

[39] Y. Hu, Z. Shi, and M. Yor (1999). Rates of convergence of diffusions with 
drifted Brownian potentials. Trans. Amer. Math. Soc. 351, 3915-3934. 

[40] N. Ikeda and H. Matsumoto (1999). Brownian motion on the hyperbolic 
plane and Selberg trace formula. J. Fund. Anal. 163, 63-110. 

[41] M. Jeanblanc, J. Pitman, and M. Yor (2002). Self-similar processes with in- 
dependent increments associated with Levy and Bessel. Stochastic Process. 
Appl. 100, 223-231. MR1919614 

[42] K. Kawazu and H. Tanaka (1993). On the maximum of a diffusion process 
in a drifted Brownian environment. In: Seminaire de Probabilites XXVII 
pp. 78-85, Lecture Notes in Math. 1557, Springer, Berlin. 

[43] K. Kawazu and H. Tanaka (1997). A diffusion process in a Brownian envi- 
ronment with drift. J. Math. Soc. Japan 49, 189-211. 

[44] V. Kac and P. Cheung (2002). Quantum Calculus. Springer, Universitext, 
New York. 

[45] H. Kesten (1973). Random difference equations and renewal theory for 
products of random matrices. Acta Math. 131, 207-248. 

[46] A. Lachal (2003). Some probability distributions in modelling DNA repli- 
cation. Ann. Appl. Probab. 13, 1207-1230. 

[47] J. W. Lamperti (1972). Semi-stable Markov processes. Z. Wahrschein- 
lichkeitstheorie verw. Gebiete 22, 205-225. 

[48] D. K. Lubensky and D. R. Nelson (2002). Single molecule statistics and 
the polynucleotide unzipping transition. Phys. Rev. E 65 03917. 

[49] H. Matsumoto and M. Yor (2000). An analogue of Pitman's 2M ~ X the- 
orem for exponential Wiener functionals. I. A time-inversion approach. 
Nagoya Math. J. 159, 125-166. 

[50] H. Matsumoto and M. Yor (2001). An analogue of Pitman's 2M - X the- 
orem for exponential Wiener functionals. II. The role of the generalized 
inverse Gaussian laws. Nagoya Math. J. 162, 65-86. 

[51] H. Matsumoto and M. Yor (2003). On Dufresne's relation between the prob- 
ability laws of exponential functionals of Brownian motions with different 
drifts. Adv. in Appl. Probab. 35, 184-206. 

[52] K. Maulik and B. Zwart (2004). Tail asymptotics for exponential fmiction- 
als of Levy processes. Stochastic Process. Appl. (to appear). 

[53] O. Mejane (2003). Equations aux differences aleatoires: des fonctionnelles 
exponentielles de Levy aux polymeres diriges en environnement aleatoire. 
These de Doctoral de I'Universite Paul Sabatier Toulouse IIL 

[54] T. Nilsen and J. Paulsen (1996). On the distribution of a randomly dis- 
counted compound Poisson process. Stochastic Process. Appl. 61, 305-310. 
[MR1386179 

[55] J. Paulsen (1993). Risk theory in a stochastic economic environment. 

Stochastic Process. Appl. 46, 327-361. 
[56] J. Paulsen and H. K. Gjessing (1997). Ruin theory with stochastic return 

on investments. Adv. in Appl. Probab. 29, 965-985. 



J. Bertoin and M. Yor / Exponential functionals 



212 



[57] M. PoUak and D. Siegmund (1985). A diffusion process and its applications 
to detecting a change in the drift of Brownian motion. Biometrika 72, 267- 
280. 

[58] V. Rivero (2005). Recurrent extensions of self-similar Markov processes and 
Cramer's condition. Bernoulli 11, 471-509. 

[59] K.-I. Sato (1991). Self-similar processes with independent increments. 
Probab. Theory Related Fields 89, 285-300. MRU 13220 

[60] K.-I. Sato (1999). Levy processes and infinitely divisible distributions. Cam- 
bridge University Press, Cambridge. 

[61] Z. Shi (2001). Sinai's walk via stochastic calculus. In: Milieux aleatoires , 
F. Comets et al. (Eds) Panoramas et Syntheses 12, pp. 53-74. 

[62] B. Simon (1998). The classical moment problem as a self adjoint finite 
difference operator. Adv. in Maths. 137, 82-203. 

[63] T. J. Stieltjes (1894-95) Recherches sur les fractions continues. Ann. Fac. 
Sci. Toulouse 8, J76-J122 and 9, A5-A47. Reprinted in : Ann. Fac. Sci. 
Toulouse Math. (6) 4 (1995). 

[64] J. Stoyanov (2000). Krein condition in probabilistic moment problems. 
Bernoulli 6, 939-949. 

[65] K. Urbanik (1992). Functionals of transient stochastic processes with inde- 
pendent increments. Studia Math. 103 (3), 299-315. 

[66] W. Vervaat (1979). On a stochastic difference equation and a representation 
of nonnegative infinitely divisible random variables. Adv. in Appl. Probab. 
11, 750-783. 

[67] D. Williams (1979). Diffusions, Markov processes, and martingales vol. 1: 

Foundations. Wiley, New York. 
[68] M. Yor (1992). Sur certaines fonctionnelles exponentielles du mouvement 

brownien reel. J. Appl. Probab. 29 , 202-208. EngHsh translation in [THj . 
[69] M. Yor (1997). Some aspects of Brownian motion. Part II: Some recent 

martingale problems. Lectures in Mathematics. ETH Ziirich, Birkhauser. 
[70] M. Yor (2001). Exponential functionals of Brownian motion and related 

processes. Springer Finance, Berlin. 



